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The electronic and magnetic properties of a neutral substitutional nickel (Ni^) impurity in dia- 
mond are studied using density functional theory in the generalized gradient approximation. The 
spin-one ground state consists of two electrons with parallel spins, one located on the nickel ion in 
the 3d^ configuration and the other distributed among the nearest-neighbor carbons. The exchange 
interaction between these spins is due to p — d hybridization and is controllable with compressive 
hydrostatic or uniaxial strain, and for sufficient strain the antiparallel spin configuration becomes 
the ground state. Hence, the Ni impurity forms a controllable two-electron exchange- coupled system 
that should be a robust qubit for solid-state quantum information processing. 



Recent advances in single-spin control of electron spins 
in quantum dots[l and spins associated with single 
dopants in semiconductors [2 suggest such systems could 
make good qubits for quantum information processing [3]. 
Spin centers based on single impurities or impurity com- 
plexes in diamond have been extensively explored as 
qubit candidates due to effective optical access and ex- 
tremely long room-temperature spin coherence times [4- 
[8]. Although a multitude of spin centers have been ob- 
served or predicted to occur in diamond, including tran- 
sition metal dopants like nickel, the spin-one nitrogen- 
vacancy (NV~) center has drawn the most attention due 
to the successful demonstration of optical spin initializa- 
tion and readout [4, 9-12 and spin manipulation via in- 
teraction with a neighboring nitrogen nuclear spin [131 [T4] 
or electron spin [15]. Ni can be introduced in diamond 
by ion implantation or during chemical vapor deposition 
[16l [It] . Characterization by optical spectroscopy, elec- 
tron paramagnetic resonance (EPR) and magnetic circu- 
lar dichroism (MCD) [T8'-'2r has revealed the existence 
of several Ni related optical centers, in particular a 1.4 eV 
optical feature which has been attributed to a Ni+ inter- 
stitial with a 3dP configuration and spin S = ^. To our 
knowledge, no experimental study of the substitutional 
Ni spin center has been reported. 

In other material systems the coherence times and fi- 
delity for quantum operations of qubits has been im- 
proved dramatically through the use of decoherence-free 
subspaces of exchange-coupled spins to form a qubit; 
one example is the use of the spin singlet state and the 
= state of a spin triplet to form a decoherence-free 
subspace for two electrons confined to two neighboring 
quantum dots[22l[22. These approaches rely on control 
of the exchange interaction between two spins, such as 
by using an electrical gate to modulate the electronic 
hopping from one quantum dot to another. Electric con- 
trol of the exchange interaction is challenging for a spin 
center in diamond, for the typical size of the electronic 
wave functions is very small, and gating technology is 
not well advanced. Strain provides an alternate mecha- 



nism for controlling a spin center, as was recently shown 
via electrically-detected magnetic resonance induced by 
strain control of the hyperfine constant of a ^^P+ donor 
in strained Si [24]. 

Here we show, using density functional theory calcu- 
lations, that the substitutional nickel impurity in dia- 
mond can be understood as an exchange-coupled system 
of two electron spins: one localized on the nickel ion and 
one delocalized on the four nearest-neighbor (NN) carbon 
atoms. The electronic configuration of Ni is unambigu- 
ously determined by p-d hybridization between the Ni 
3d and NN carbon 2p levels. Although the ground state 
at ambient pressure for the neutral nickel impurity has 
been predicted to be a spin-one center [25] , the spin-zero 
state is nearly degenerate, and can be made degenerate 
through the application of reasonable compressive hy- 
drostatic or uniaxial strain. To reduce to an effective 
two-state decoherence free subspace, similar to that im- 
plemented for quantum dots [23 , the = ±1 triplet 
states (T±) could be split off by a magnetic field. For the 
nickel ion, strain modulation could be used to manipu- 
late the energy splitting between the singlet (S) and the 
remaining Sz = triplet state (Tq), instead of the electro- 
static gating used for double quantum dots [23 . As the 
environments differ for the two electron spins that com- 
bine to form the S or Tq states, the slightly different g 
factors we find for these two spins provide an orthogonal 
axis of control in the effective two-state S/Tq subspace 
of the nickel ion. 

The calculations were performed with the scalar 
relativistic version of the full potential local orbital 
FPL09. 00-33 code [26] using the (spin-polarized) gen- 
eralized gradient approximation ((S)GGA) with the 
parametrization of Perdew, Burke and Ernzerhof [27J. 
The convergence of the results with respect to k-space in- 
tegrals was carefully checked, and 8x8x8 = 512 k-points 
were sufficient. We used a 64- atom NiCea supercell corre- 
sponding to a 2 X 2 X 2 multiple of the (non-fundamental) 
cubic unit cell of diamond with Ni substituted for one car- 
bon. The supercell size was fixed to correspond to the ex- 
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perimental lattice constant of diamond ao = 3.5668A[28^, 
and all the atomic positions within the supercell were 
allowed to relax with a precision of 1 meV/A. We con- 
sidered two possible symmetries for the Ni impurity, 
with four identical NN carbons, and Csy in which a trig- 
onal distortion along the [111] axis is allowed. We found 
that the relaxation primarily involves the NN carbons 
which moved by 14.2% while further neighbors moved by 
< 1.1%. The Csy relaxation along [111] was less than 
0.05% and so we considered only in subsequent calcu- 
lations. The electron occupation number of Ni was 26.9 
and the total magnetization of the ground state was 2.0 
/iB, with approximately 0.8 fiB localized on the Ni, and 
the rest distributed among the NN carbons. 
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FIG. 1. Total and Ni partial SGGA DOS of a NiCes supercell. 

Fig. [T] shows the total and Ni partial SGGA density 
of states (DOS) of the NiCes supercell. The calculated 
direct band gap of 4.9 eV and indirect gap of 4.8 eV are 
both smaller than the experimental gap Eq"^ = 5.47 eV, 
due to the gap underestimation of the GGA functional 
[29] . The Ni 3d levels are split by the crystal field into 
doubly degenerate Cg states and a higher energy triplet 
of t2g levels. The spin up and spin down Cg levels are lo- 
calized on the Ni site, giving rise to peaks approximately 
1 eV below the valence band maximum (VBM) and sep- 
arated by a Hund exchange splitting = 0.4 eV. The 
t2g levels strongly hybridize with the valence band (pre- 
dominantly 2p character), forming a partially filled and 
a totally empty bound state at 2.5 and 3.0 eV above the 
VBM respectively, corresponding to the spin up and spin 
down antibonding levels as described below. We note 
the presence of a ferromagnetic (FM) spin splitting of 
the valence band corresponding to the strong hybridiza- 
tion limit similar to the case of Mn^+ in ZnO or Mn^+ 
in GaN [SOlET]. 

The separate electronic configurations of the Ni and 
the neighboring carbon atoms clarify the electronic na- 



ture of the impurity. Ni is in the positively charged 
[Ar]4s^3(i^ configuration Ni+ with a spin Si = ^. 
Nickel's positive charge induces a single electron in the 
surrounding diamond that is distributed among the dan- 
gling bonds of the NN carbons. This addition of the in- 
duced electron to the four dangling bonds, each with 2^2^ 
character, results in a 2s^2p^ configuration with S2 = \- 
The levels of the Ni+ and NN carbons hybridize, as 
shown in Fig. [2| The 2p derived ^2 defect level hy- 
bridizes with the Ni t2g levels to form bonding (t^ ) and 
antibonding (t^^) levels. This p-d hybridization picture 
identifies the origin of the calculated SGGA DOS (Fig. 
ni) as due to a FM interaction between and 82- 




FIG. 2. Hybridization between the Ni Cg and t2g Sd levels 
and the 2p-derived nearest-neighbor carbon dangling bond 
with a ferromagnetic alignment providing a total spin 5't = 1 
(triplet state T). 
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FIG. 3. Heisenberg exchange coupling between the two spins 
5*1 and S2 calculated as a function of hydrostatic strain. 

Strain dramatically modifies the exchange coupling be- 
tween the spins of the Ni+ and the surrounding carbons. 
The energy difference between the FM and AFM ar- 
rangements of the two spins was equated to the total 
energy difference between the SGGA and GGA calcula- 
tions = {EsGGA - Egga)/2 = -{J/2)St{St + 1), 
with St = Si^ S2 and the Hamiltonian H = — 2JSi • S2. 
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Hydrostatic strain was included by changing the size of 
the supercell, within which the atomic positions were al- 
lowed to relax. For [001] uniaxial strain, we fixed the 
atomic positions to their calculated positions for the un- 
strained system and rescaled the supercell anisotropically 
according to the Poisson ratio [32] . 

Fig. [3] shows the exchange coupling as a function of 
hydrostatic strain, where J > and J < correspond to 
FM and AFM coupling respectively. The ground state 
has an AFM spin alignment for en < —7 %. As a check, 
we used an initial configuration of 1 /ie on the Ni and 
-1 fiB on the NN carbons to locate the low-spin (non- 
ground-state) AFM solution with the SGGA functional 
for en = 0. The total energy calculated for this solu- 
tion is within 3 meV of the nonmagnetic GGA solution, 
which is within our estimated error. The SGGA AFM 
ground state for en = —'^0 % is similarly degenerate with 
the GGA result, and the calculated total magnetization 
is nearly zero (~ 0.1 jj^b)- Under [001] uniaxial strain, 
the ground state becomes AFM for a compressive strain 
6[ooi] < —19 %• The existence of a pressure induced tran- 
sition between a high-spin and low-spin state is also ob- 
served in transition- metal-ion compounds, where it was 
shown theoretically to be due to a competition between 
localization (induced by the Hund's exchange coupling 
which favors the high-spin state) and a tendency towards 
delocalization induced by the crystal field (which favors 
the low-spin state) [33] . 
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FIG. 4. Total and Ni partial SGGA DOS of < 
with a compressive hydrostatic strain en = 



NiCes supercell 
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Fig. |4] presents the total and Ni partial SGGA DOS of 
the AFM ground state with en = —'^0 %. We found that 
the band gap increases to 5.8 eV. The spin up and spin 
down eg levels are now spin degenerate (J^ = 0), and 
there is a hybridization between the Ni t2g and the NN 
C 2p levels, leading to the formation of bonding and an- 
tibonding levels. The spin-degenerate antibonding levels 
are partially filled at the Fermi level and correspond to 
a bound state in the gap 3.4 eV above the VBM. Fig. 



|5] gives a schematic representation of the p-d hybridiza- 
tion model under strain. The crystal field splits the Ni 
3d level into lower Cg and upper t2g levels which are now 
spin degenerate (J^ =0)- To explain the non-magnetic 
DOS presented in Fig. [4| we need to take into account 
an AFM interaction between the two spins and 5*2 . 
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FIG. 5. Electronic configuration corresponding to an AFM 
alignment of Si and S2 , providing a total spin 5't = (singlet 
state S). 

The Ni^ system can be described in the singlet S 
(AFM, St = 0) - triplet T (FM, St = 1) basis by the 
Hamiltonian = -2J{{e})Si-S2-^/J^B{9iSiz-^g2S2z)Bz- 
With the T± states split off by an applied magnetic field 
Bz^ we have a two- level S/Tq system (qubit). In the S/Tq 
subspace, H can be written as follows: 



/ £s{{e}) Ag^ 
\ AG^ 



(1) 



with £s{{e}) = 2J({e}) the strain-controlled exchange 
splitting {e.g. Fig. |3|, AQ^ = {jiiBBzAg)/4 and the basis 
states \S) = ^(1 tNilc) - | iNitc)), |To) = ^(| tNilc) 
+ 1 iNitc))- Rapid dynamical spin exchange occurs be- 
tween the Ni and C sites, with a speed that can be es- 
timated as ~ 1 fs from the ~ 3 eV splitting between 
bonding and antibonding t2g states in Fig. |2] 

The singlet-triplet mixing AQ^ depends on the differ- 
ence in the Lande factors Ag = gi — g2 between the spin 
Si localized on the Ni and the spin ^2 located on the NN 
carbon atoms. Applying Hund's rule to the Ni spin we 
obtain a gi factor of 1.33, whereas the spin on the NN car- 
bons is delocalized and thus should be well described by a 
free-electron ^2 factor of 2. We obtain a finite difference 
Ag ~ —0.66. The Ni^ qubit, a two exchange-coupled 
spin-^ system, can therefore be used to develop an uni- 
versal quantum computer (QC) architecture [22]. The 
spins Si and 5*2 form nearly decoherence free subspaces 
(DFS) due to the long spin lifetime in diamond. Based 
on these DFSs, a solid-state QC architecture based on 
Ni^ qubit should be immune from collective decoherence 
mechanism. 

Generation of sufficient strain to make the FM (Tq) 
and AFM (S) states degenerate would require a stress of 
^ 80 GPa, which is well below the > 300 GPa maximum 
of a diamond anvil cell [34]; such strains and stresses are 
also typical in nanoscale pseudomorphic structures, such 
as InAs/InSb heterostructure nanowires [35 . Once the S 
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and To states are degenerate, high-speed manipulation of 
the two-state system with modulated strain can be done 
with strain modulation amplitudes orders of magnitude 
less. To achieve an energy splitting corresponding to a 
~ 1 GHz precession frequency, a strain modulation of 
only 10~^ would be required, which is three orders of 
magnitude smaller than that used in Ref. High- 
Q nanomechanical diamond resonators operating above 
1 GHz have already been demonstrated [36]. 
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FIG. 6. Ni° qubit (top) : singlet S (AFM) - triplet T (FM) 
groundstate of the spins 5*1 and 5*2 (black arrows) depending 
on strain, qubit initialization and potential V(R) for the two 
electronic spins as a function of their separation R (bottom). 



Fig. [g] shows (top) schematically the location and ori- 
entation of the spins in the Ni^ qubit system. The Ni 
of spin = ^ is at the center of the tetrahedra formed 
by the four NN carbon atoms with a spin 6*2 = ±^ dis- 
tributed over the four NN carbon dangling bonds, de- 
scribing a pseudo-orbit around the Ni atom. Whereas 
in the simplest picture of quantum-dot gate control it 
is the barrier between the spins which is lowered, here 
the potential minima are brought closer together and the 
exchange splitting £s is controlled by the application of 
strain (Eq. [T]). 

The qubit can be initialized in the singlet state (e = e^, 
St = 0) by the application of strain (Fig. [6|, or in the T+ 
triplet state (e = ex, ^S't = 1) in a static magnetic field. 
Figure [l] also reveals the potential of spin-selective opti- 
cal processes to initialize the spin. Spin-orbit interaction 
within the partially-occupied spin-polarized t2g manifold 
of the Ni ion, not shown in the diagrams so far, will split 
the three states according to the orbital angular momen- 
tum projection of the t^g state (+1, 0, —1) parallel to the 
spin, with the state with spin and orbit parallel at the 
highest energy (and thus unoccupied). Optical pump- 



ing with circularly-polarized light at about 3 eV, corre- 
sponding to a transition from the valence maximum to 
the lowest-energy unoccupied minority-spin t2g feature 
in the density of states, will create minority spins in the 
t2g states if minority spin direction and light polarization 
are antiparallel, but not if parallel. Eventually, non-spin- 
selective recombination combined with spin-selective op- 
tical generation of carriers will drive the ground state 
spin = 1 to point antiparallel to the polarization of 
the optical pump. From the state a microwave pulse 
can coherently manipulate the spin into the Tq state. 
Once initialized, the qubit can be manipulated within the 
/S/To DFS by strain modulation and resonant microwave 
radiation [9 through Eq. ([T]). 

Readout of the quantum state of the Ni^ qubit could 
be performed by electron paramagnetic resonance, or 
through spin-selective optical measurements. Because 
of the spin-orbit interaction described above, an opti- 
cal probe with photon energy just below the valence to 
unoccupied t^g state transition will experience Faraday 
rotation that depends on the orientation of the Ni spin, 
which allows a measurement of the state of the system in 
the I tNilc), I iNitc) basis. 

In conclusion, we have performed ah initio calculations 
of the ground state of Ni^ in diamond. The individual 
spins associated with this dopant, that of the Ni and 
of the nearest-neighbor C's, should each be limited in 
coherence time by similar processes in diamond that limit 
the coherence time of the NV~ spin center, especially the 
presence of nuclear spins (^^Ni nuclear spins have similar 
rarity to ^^C nuclear spins). Thus, due to the use of 
a decoherence-free subspace, a qubit architecture based 
on Ni^ may be more immune from decoherence than a 
double-quantum-dot-based spin qubit or an NV~ spin 
center in diamond. 
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